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Description 



METHOD AND APPARATUS FOR MEASUREMENT AND CORRECTION OF REFRACTIVE POWER 
^DISTRIBUTION DATA : : : — : ~ " 



Technical Field 

The present invention relates to a correction data 
10 measurement method for measuring appropriate correction data, 
a measuring apparatus, and a recording medium recording a 
measurement program. 

Background Art 

15 

Conventionally, as a technique for measuring ocular 
correction data, measurement of S (Sphere), C (Cylinder) and 
A (axis) by a ref ractometer has been carried out. Besides, 
recently, an eye characteristic measuring apparatus capable of 

20 measuring higher order aberrations has also been developed, and 
not only S, C and A on a line like, for example, a ring of <j>3 
mm as in a ref ractometer , but also S, C and A on a plane when 
a pupil diameter is made various sizes can be calculated from 
lower order aberrations. By the eye characteristic measuring 

25 apparatus like this, especially after a refraction correcting 
surgical operation or in an eye disease, values closer to 
prescription values of eyeglasses or contact lenses than the 
ref ractometer can be calculated (for example, see Japanese 
Patent Application No. 2001-119145, No. 2001-120002, No. 
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2001-119086, No. 2000-318534, etc.). 

However, there ±s still a case where a difference occurs 
between the objective calculation result of a conventional eye 
characteristic measuring apparatus and the prescription value 
5 of eyeglasses, contact lenses, lenses or the like, and there 
has been a case where it is insufficient as the evaluation of 
S, C and A. 

Disclosure of the Invention 

10 

Accordingly, in the present invention, in a case where 
higher order aberration are included in a result of a 
measurement obtained by an eye characteristic measuring 
apparatus capable of measuring higher order aberrations , lower 

.15 order aberrations corresponding to the time of an objective 
complete correction are not made corrective correction data, 
optical performance is evaluated by, for example, a Strehl 
ratio or a phase shift, such lower order aberration quantities 
that the Strehl ratio are high and/or the phase shift is 

20 decreased is calculated, and corrective correction data of S, 
C, A and the like at that time is obtained, so that more optimum 
correction data close to a prescription value of eyeglasses, 
contact lenses or the like is obtained, which is an object of 
the invention. 

25 According to first solving means of the invention, a 

correction data measurement method comprises 

a first step at which an arithmetic part receives 
measurement data indicating a refractive power distribution of 
a subject eye and obtains lower order aberrations and higher 

2 
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order aberrations on the basis of the measurement data, 

a second step at which the arithmetic part judges whether 
the higher order aberrations have a specified value or higher , 
and 

5 a third step at which the arithmetic part changes, in a 

case where the higher order aberrations have the specified 
value or higher, lower order aberration quantities 
corresponding to the higher order aberrations having the 
specified value or higher and obtains appropriate correction 

10 data suitable for the subject eye. 

Besides, there are provided a correction data 
measurement program for causing a computer to execute the above 
respective steps, and a computer readable recording medium 
recording the correction data measurement program. 

15 According to second solving means of the invention, a 

correction data measuring apparatus comprises an arithmetic 
part for obtaining an optical characteristic of a subject eye 
by performing a Zernike analysis on the basis of inclination 
angles of light fluxes obtained by a first light receiving part, 

20 wherein 

the arithmetic part includes 

first means for receiving measurement data indicating a 
refractive power distribution of the subject eye and obtaining 
lower order aberrations and higher order aberrations on the 
25 basis of the measurement data, 

second means for judging whether the higher order 
aberrations have a specified value or higher, and 

third means for changing, in a case where the higher order 
aberrations have the specified value or higher, lower order 
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aberration quantities corresponding to the higher order 
aberrations having the specified value or higher and obtaining 
appropriate correction data suitable for the subject eye. 

According to third solving means of the invention, a 
5 correction data measurement method comprises steps of: 

receiving measurement data indicating a refractive power 
distribution of a subject eye and obtaining lower order 
aberrations and higher order aberrations on the basis of the 
measurement data, 
10 forming a relational expression of a Strehl ratio and a 

phase shift (PTF) from the obtained lower order aberrations and 
higher order aberrations, and 

changing the lower order aberration to obtain a condition 
under which the Strehl ratio becomes maximum and the phase shift 
15 (PTF) becomes substantially zero, and making lower order 
aberration quantities at that time a correction value. 

Brief Description of the Drawings 

20 Fig, 1 is a structural view of an optical system 100 of 

an eye optical characteristic measuring apparatus. 

Fig. 2 is a structural view of an electrical system 200 
of the eye optical characteristic measuring apparatus. 

Fig. 3 is an explanatory view of a Landolt's ring. 
25 Fig. 4 is a flowchart of an eye characteristic 

measurement . 

Fig. 5 is a flowchart of a correction image simulation 
of step S107. 

Fig. 6 is a flowchart of a first example of a best image 
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condition calculation. 

Fig* 7 is a flowchart of a second example of a best image 
condition calculation. 

Fig. 8 is a view showing a display example concerning best 
5 image display - Strehl ratio optimization (pupil diameter of 
, 4 mm) . - . - - . - 

Fig. 9 is a view showing a display example concerning best 
image display -r PTF optimization (pupil diameter of 4 mm) . 

Fig. 10 is a view showing a display example concerning 
10 the comparison (pupil diameter of 4 mm) of corrections added 
with a lower order correction and a higher order correction. 

Fig. 11 is an explanatory view of an example of 
prescription data (pupil diameter of 4 mm) for eyeglasses or 
contact lenses. 

15 Fig. 12 is an explanatory view of an example of data 

(pupil diameter of 4 mm) for a refraction correcting surgical 
operation. 

Fig. 13 is a view of Zernike coefficients of (r, t) 
coordinates . 

20 Fig. 14 is a view of Zernike coefficients of (x, y) 

coordinates . 

Best Mode for Carrying Out the Invention 
25 1. Eye optical characteristic measuring apparatus 



Fig. 1 is a structural view of an optical system 100 of 
an eye optical characteristic measuring apparatus. 

The optical system 100 of the eye optical characteristic 
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measuring apparatus is an apparatus for measuring an optical 
characteristic of an eye 60 to be measured as an object, and 
includes a first illuminating optical system 10, a first light 
receiving optical system 20, a second light receiving optical 
5 system 30, a common optical system 40, an adjusting optical 
system 50 , a second illuminating optical system 70 , and a second 
light sending optical system 80. Incidentally , with respect to 
the eye 60 to be measured, a retina 61 and a cornea 62 are shown 
in the drawing. 

10 The fiyst illuminating optical system 10 includes, for 

example, a first light source part 11 for emitting a light flux 
of a first wavelength, and a condensing lens 12, and is for 
illuminating a minute area on the retina (retina) 61 of the eye 
60 to be measured Vith the light flux from the first light source 

15 part 11 so that its illumination condition can be suitably set. 
Incidentally, here, as an example, the first wavelength of the 
illuminating light flux emitted from the first light source 
part 11 is a wavelength (for example, 780 nm) of an infrared 
range . 

20 Besides, it is desirable that the first light source part 

11 has a high spatial coherence and a low temporal coherence. 
Here, the first light source part 11 is, for example, a super 
luminescence diode (SLD) , and a point light source having high 
luminescence can be obtained. Incidentally, the first light 

25 source part 11 is not limited to the SLD, and for example, a 
laser having a high spatial coherence and a high temporal 
coherence can also be used by inserting a rotation diffused 
plate or the like to suitably lower the temporal coherence. 
Further, an LED having a low spatial coherence and a low 
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temporal coherence can also be used, if light quantity is 
sufficient, by inserting, for example, a pinhole or the like 
at a position of a light source in "an optical path. 

The first light receiving optical system 20 includes, for 
5 example, a collimator lens 21, a Hartmann plate 22 as a 
conversion member for converting a part of a light flux (first 
light flux) reflected and returned from the retina 61 of the 
eye 60 to be measured into at least 17 beams, and a first light 
receiving part 23 for receiving the plural beams converted by 
10 the Hartmann plate 22, and is for guiding the first light flux 
to the first light receiving part 23. Besides, here, a CCD with 
little readout noise is adopted for the first light receiving 
part 23, and as the CCD, a suitable type of CCD, for example, 
a general low noise type of CCD, a cooling CCD of 1000 * 1000 
15 elements for measurement, or the like is applicable. 

The second illuminating optical system 70 includes a 
second light source 72 and a Placido f s disk 71. Incidentally, 
the second light source 72 can be omitted. The Placido's disk 
(PLACIDO f S DISK) 71 is for projecting an index of a pattern 
20 composed of plural co-axial rings. Incidentally, the index of 
the pattern composed of the plural co-axial rings is an example 
of an index of a specified pattern, and a different suitable 
pattern can be used. Then, after an alignment adjustment 
described later is completed, the index of the pattern composed 
25 of the plural co- axial rings can be projected. 

The second light sending optical system 80 is for mainly 
performing, for example, the alignment adjustment described 
later, and measurement and adjustment of a coordinate origin 
and a coordinate axis, and includes a second light source part 
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31 for emitting a light flux of a second wavelength, a 
condensing lens 32, and a beam splitter 33. 

The second light receiving optical system 30 includes a 
condensing lens 34 and a second light receiving part 35 . The 
5 second light receiving optical system 30 guides a light flux 
(second light flux), which is originated from the pattern of 
the Placido's disk 71 illuminated from the second illuminating 
optical system 70 and is reflected and returned from the 
anterior eye part or the cornea 62 of the eye 60 to be measured, 
10 to the second light receiving part 35. Besides, it can also 
guide a light flux, which is emitted from the second light 
source part 31 and is reflected and returned from the cornea 
62 of the eye 60 to be measured, to the second light receiving 
part 35. Incidentally, as the second wavelength of the light 
15 flux emitted from the second light source part 31, for example, 
a wavelength different from" the first wavelength (here, 780 nm) 
and longer (for example, 940 nm) than that can be selected. 

The common optical system 40 is disposed on an optical 
axis of the light flux emitted from the first illuminating 
20 optical system 10, can be included in the first and the second 
illuminating optical systems 10 and 70, the first and the second 
light receiving optical systems 20 and 30, the second light 
sending optical system 80 and the like in common, and includes, 
for example, an afocal lens 42, beam splitters 43 and 45, and 
25 a condensing lens 44. The beam splitter 43 is formed of such 
a mirror (for example, a dichroic mirror) that the wavelength 
of the second light source part 31 is sent (reflected) to the 
eye 60 to be measured, and the second light flux reflected and 
returned from the retina 61 of the eye 60 to be measured is 

8 
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reflected, and on the other hand, the wavelength of the first 
light source part 11 is transmitted. The beam splitter 45 is, 
formed of such a mirror (for .example, a polarization beam ^ 
splitter) that the light flux of the first llghtf source part 
5 11 is sent (reflected) to the eye 60 to be measured, and the 
first light flux reflected and returned from the retina 61 bf 
the eye 60 to be measured is transmitted. By the beam splitters 
43 and 45, the first and the second light fluxes do not mutually 
enter the other optical systems to generate noise. 

10 The adjusting optical system 50 is for 'mainly performing, 

for example, a working distance adjustment described later, 
includes a third light source part 51, a fourth light source 
part 55, condensing lenses 52 and 53, and a third light 
receiving part 54, and is for mainly performing the working 

15 distance adjustment. 

The third illuminating optical system 90 includes an 
optical path for projection of an index for causing, for example, 
fixation of the subject eye or fogging, and includes a fifth 
light source part (for example, a lamp) 91, a fixed index 92 

20 and a relay lens 93. The fixed index 92 can be irradiated to 
the retina 61 by the light flux from the fifth light source part 
91, and the subject eye 60 is made to observe its image. The 
fixed index 92 and the retina 61 are put in a conjugated relation 
by the third illuminating optical system 90. 

25 A sixth driving part 915 is for moving, for example, the 

fixed index 92 of the third illuminating optical system 90, and 
outputs a signal (15) to not-shown suitable movement means and 
drives this movement means. By this, the sixth driving part 915 
can perform the movement and adjustment of the fixed index 92 
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of the third Illuminating optical system 90. 

Next, the alignment adjustment will be described* The 
alignment adjustment is mainly carried but by the second light 
receiving optical system 30 and the* second light sending 
5 optical system 80. 

First, the light flux from the second light source part 
31 illuminates the eye 60 to be measured as the object with the 
substantially parallel light flux through the condensing lens 
32, the beam splitters 33 and 43, and the afocal lens 42. The 

10 reflected light flux reflected by the cornea 62 of the eye 60 
to be measured is emitted as a divergent light flux such as is 
emitted from a point at the half of the radius of curvature of 
the cornea 62. The divergence light flux is received as a spot 
image by the second light receiving part 35 through the afocal 

15 lens 42, the beam splitters 43 and 33, and the condensing lens 
34. . 

Here, in the case where the spot image on the second light 
receiving part 35 is outside the optical axis, the main body 
of the eye optical characteristic measuring apparatus is moved 

20 and adjusted vertically and horizontally, and the spot image 
is made to coincide with the optical axis . As stated above , when 
the spot image coincides with the optical axis, the alignment 
adjustment is completed. Incidentally, with respect to the 
alignment adjustment, the cornea 62 of the eye 60 to be measured 

25 is illuminated by the third light source part 51, and an image 
of the eye 60 to be measured obtained by this illumination is 
formed on the second light receiving part 35, and accordingly,, 
this image may be used to make the pupil center coincide with 
the optical axis. 

10 
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Next, the working distance adjustment will be described. 
The working distance adjustment is mainly carried out by the 
adjusting optical system 50. 

First the working distance adjustment is carried out by, 
5 for example, irradiating the eye 60 to be measured with a 
parallel light flux emitted from the fourth light source part 
55 and close to the optical axis, and by receiving the light 
reflected from the eye 60 to be measured through the condensing 
lenses 52 and 53 by the third light receiving part 54. Besides, 

10 in the case where the eye 60 to be measured is in a suitable 
working distance, a spot image from the fourth light source part 
55 is formed on the optical axis of the third light receiving 
part 54. On the other hand, in the case where the eye 60 to be 
measured goes out of the suitable working distance, the spot 

15 image from the fourth light source part 55 is formed above or 
below the optical axis of the third light receiving part 54 . 
Incidentally, since the third light receiving part 54 has only 
to be capable of detecting a change of a light flux position 
on the plane containing the fourth light source part 55, the 

20 optical axis and the third light receiving part 54 , for example, 
a one -dimensional CCD arranged on this plane, a position 
sensing device (PSD) or the like is applicable. 

Next, a positional relation between the first 
illuminating optical system 10 and the first light receiving 

25 optical system 20 will be described. 

The beam splitter 45 is inserted in the first light 
receiving optical system 20, and by this beam splitter 45, the 
light from the first illuminating optical system 10 is sent to 
the eye 60 to be measured, and the reflected light from the eye 

11 
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60 to be measured is transmitted. The first light receiving part 
. 23 included in the first light receiving optical system 20 
receives the light transmitted through the Hartmann plate 22 
as the conversion member and generates a received light signal . 
5 Besides, the first light source part 11 and the retina 

v 61 of the' eye 60 to be measured form a conjugated relation. The 
retina 61 of the eye 60 to be measured and the first light 
receiving part 23 are conjugate. Besides , the Hartmann plate 
22 and the pupil of the eye 60 to be measured form a conjugated 
10 relation. Further, the first light receiving optical system 20 
forms a substantially conjugated relation with respect to the 
cornea 62 as the anterior eye part of the eye 60 to be measured,* 
the pupil, and the Hartmann plate 22. That is, the front focal 
point of the afocal lens 42 is substantially coincident with 
15 the cornea 62 as the anterior eye part of the eye 60 to be 
measured and th&~ pupil. 

Besides, the first illuminating optical system 10 and the 
first light receiving optical system 20 are moved together so 
that a signal peak according to the reflected light at the light 
20 receiving part 23 becomes maximum on the condition that the 
light flux from the first light source part 11 is reflected at 
a point on which it is condensed. Specifically, the first 
illuminating optical system 10 and the first light receiving 
optical system 20 are moved in a direction in which the signal 
25 peak at the first light receiving part 23 becomes large, and 
are stopped at a position where the signal peak becomes maximum. 
By this, the light flux from the first light source part 11 is 
condensed on the eye 60 to be measured. 

Besides, the lens 12 converts a diffused light of the 

12 
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light source 11 into a parallel light. A diaphragm 14 is 
positioned at an optically conjugated position with respect to 
the pupil of the eye or' the Hartmann^ plate 22. The diaphragm 

14 has a diamfeter smaller than an effective range of the 
5 Hartmann plate 22, and the so-called single path aberration 

measurement (method 4 in which aberrations of an eye have an 
influence on only the light receiving side) is established. In 
order to satisfy the above, the lens 13 is disposed such that 
the retina conjugated point of the real light beam coincides 

10 with the front focal position, and further, in order to satisfy 
the conjugated relation between the lens and the pupil of the 
eye, it is disposed such that the rear focal position coincides 
with the diaphragm 14. 

Besides, after a light beam 15 comes to have a light path 

15 common to a light beam 24 by the beam splitter 45, it travels 
in the same way as the light beam 24 paraxially. However, in 
the single path measurement, the diameters of the light beams 
are different from each other, and the beam diameter of the 
light beam 15 is set to be rather small as compared with the 

20 light beam 24. Specifically, the beam diameter of the light beam 

15 is, for example, about 1 mm at the pupil position of the eye, 
and the beam diameter of the light beam 24 can be about 7 mm 
(incidentally, in the drawing, the light beam 15 from the beam 
splitter 45 to the retina 61 is omitted) . 

25 Next, the Hartmann plate 22 as the conversion member will 

be described. 

The Hartmann plate 22 included in the first light 
receiving optical system 20 is a wavefront conversion member 
for converting a reflected light flux into plural beams. Here, 

13 
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plural micro -Fresnel lenses disposed on a plane orthogonal to 
the optical axis apply in the Hartmann plate 22. Besides , in 
general, with respect to the measurement object part (the eye 
60 to be measured) > in order to measure a sphere of the eye 60 
5 to be measured, third- order astigmatism aberrations, and other 
higher order aberrations , it is necessary to - perform the 
measurement with at least 17 beams through the eye 60 to be 
measured . 

Besides, the micro-Fresnel lens is an optical element, 
10 and includes, for example, a ring with a height pitch for each 
wavelength, and a blade optimized for emission parallel to a 
condensing point. The micro-Fresnel lens here is subjected to, 
for example, 8 -level optical path length variation employing 
g. semiconductor fine working technique, and achieves a high 
15 condensing efficiency (for example, 98 %). 

Besides, the reflected light from the retina 61 of the 
eye 60 to be measured passes through the af ocal lens 42 and the 
collimate lens 21 and is condensed on the first light receiving 
part 23 through the Hartmann plate 22. Accordingly, the 
20 Hartmann plate 22 includes a wavefront conversion member for 
converting the reflected light flux into at least 17 beams. 

Fig. 2 is a block diagram showing an electrical system 
200 of the eye optical characteristic measuring apparatus. The 
25 electrical system 200 of the eye optical characteristic 
measuring apparatus includes, for example, an arithmetic part 
210, a control part 220, a display part 230, a memory 240, a 
first driving part 250, and a second driving part 260. 

The arithmetic part 210 receives a received light signal 

14 
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(4) obtained from the first light receiving part 23 , a received 
light signal (7) obtained from the second light receiving part 
35/ and a received light signal (10) obtained from the third 
light receiving part 54 , and performs an arithmetical operation 
5 on the origin of coordinates, a coordinate axis, movement of 
coordinates, rotation, ocular aberrations, corneal 
aberrations, Zernike coefficients , aberration coefficients , a 
Strehl ratio, a white light MTF, a Landolt's ring pattern and 
the like. Besides, signals corresponding to such calculation 
10 results are outputted to the control part 220 for performing 
the whole control of an electric driving system, the display 
part 230, and the memory 240, respectively. Incidentally, the 
details of the arithmetic part 210 will be described later. 

The control part 220 controls lighting and extinction of 
15 the first light source part 11 on the basis of the control signal 
from the arithmetic part 210, or controls the first driving part 
250 and the second driving part 260. For example, on the basis 
of the signals corresponding to the operation results in the 
arithmetic part 210, the control part outputs a signal (1) to 
20 the first light source part 11, outputs a signal (5) to the 
Placido's disk 71, outputs a signal (6) to the second light 
source part 31, outputs a signal (8) to the third light source 
part 51, outputs a signal (9) to the fourth light source part 
55, outputs a signal (11) to the fifth light source part 91, 
25 and outputs signals to the first driving part 250 and the second 
driving part 260. 

The first driving part 250 is for moving the whole first 
illuminating optical system 10 in the optical axis direction 
on the basis of, for example, the received light signal (4) 

15 
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inputted to the arithmetic part 210 from the first light 

receiving part 23, and outputs a signal (2) to a not-shown 

suitable lens movfettfent means and drives the lens % movement means . 

By this, the first driving part 250 can perform the movement 

5 and adjustment of the first illuminating optical system 10, 

The second driving part 260 is for moving the whole first 

light receiving optical system 20 in the optical axis direction 

on the basis of, for example, the received light signal (4) 

inputted to the arithmetic part 210 from the first light 

10 receiving part 23, and outputs a signal (3) to a not -shown 

♦ 

suitable lens movement means, and drives the lens movement 

means. By this, the second driving part 260 can perform the 

# 

movement and adjustment of the first light receiving optical 
system 20. 
15 

2. Zernike analysis 

Next, a Zernike analysis will be described. A generally 
known method of calculating Zernike coefficients Ci 23 " 1 from 

20 Zernike polynomials will be described. The Zernike 
coefficients Ci 23 " 1 are important parameters for grasping the 
optical characteristic of the subject eye 60 on the basis of 
inclination angles of the light fluxes obtained by the first 
light receiving part 23 through the Hartmann plate 22. 

25 Wavefront aberrations W(X, Y) of the subject eye 60 are 

expressed using the Zernike coefficients Ci 23 " 1 and the Zernike 
polynomials Zt 23 " 1 by the following expression. 
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n i 



V-> Z fJ-'(x,Y) 



i=0 7=0 



Where, (X, Y) 



denotes vertical and horizontal 



coordinates of the Hartmann plate 22. 

Besides, with respect to the wavefront aberrations W(X, 
5 Y) , when the horizontal and vertical coordinates of the first 
light receiving part 23 are denoted by (x, y) , a distance 
between the Hartmann plate 22 and the first light receiving part 
23 is denoted by f , and a movement distance of a point image 
received by the first light receiving part 23 is denoted by (Ax, 
10 Ay), the following expression is established. 



Where, the Zernike polynomials Zi. 23 " 1 are expressed by the 
following numerical expressions. Specifically, Fig. 13 is a 
15 view of the Zernike coefficients of (r, t ) coordinates, and Fig. 
14 is a view of the Zernike coefficients of (x, y) coordinates. 



dW(X,Y) 
dX 



A* 
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m>0 sin 



m^0 cos 
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Incidentally, with respect to the Zernike coefficients 
d 23 " 1 , specif ic values can be obtained by minimizing the squared 
error expressed by the following numerical expression. 



data number 



'f ay(*,.r f ) Ax,]* \ dW{X l9 Y,) Ay,' 
{ dX f) \ BY •/ J 



Where, W(X # Y) : wavefront aberrations, (X, Y): Hartmann 
plate coordinates, (Ax, Ay): a movement distance of a point 
image received by the first light receiving part 23, f: a 
distance between the Hartmann plate 22 and the first light 
10 receiving part 23 r m: the number of data. 

The arithmetic part 210 calculates the Zernike 
coefficients Ci 23-1 , and uses this to obtain eye optical 
characteristics such as spherical aberrations, coma 
aberrations, and astigmatism aberrations. 

15 

3. Landolt's ring 



Fig. 3 is an explanatory view of a Landolt's ring. 
Hereinafter, preparation of data of a luminous 
20 distribution function Land(x, y) of the Landolt's ring will be 
described. 

The Landolt's ring is expressed by the reciprocal of a 
recognizable minimum visual angle, and the ability to be 

18 
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capable of recognizing a visual angle of one minute is called 
visual acuity of 20/20. For example , if the recognizable 
minimum visual angle is 2 minutes, the visual acuity is defined 
as 20/40 , and if 10 minutes, the visual acuity is defined as 
5 20/200. In general, the Landolt's ring uses, as an index, a ring 
in which a gap being 1/5 of the size 6f the outside ring is 
provided as shown in the drawing. 

When the visual acuity is V, the size d of the Landolt's 
ring projected on the retina is calculated by 

( 1 1^ 
<*=5x2-tftan x- 

UO-F 2) 

10 

(R: a distance between a pupil and an image point (retina)) 
On the basis of this expression and the definition of the 
Landolt's ring, a black portion of the Landolt's ring is made 
0 , a white portion thereof is made 1 , and the luminous 
15 distribution function Land(x, y) of the Landolt's ring is 
prepared. The data of the prepared luminous distribution 
function Land(x, y) is stored in the memory 240, is read out 
by the arithmetic part 210, and is set correspondingly to 
predetermined visual acuity. 

20 

4. Correction data measurement method 

Fig. 4 is a flowchart of an eye characteristic 
measurement . 

25 First, the eye characteristic measuring apparatus makes 

an alignment of X, Y and Z axes of the pupil position of the 
subject eye (S101). Next, the measuring apparatus moves the 

19 
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origin of a movable part (S103) . For example, the Hartmann plate , 
the Placido's disk or the like is matched to zero diopter. The 
arithmetic part 210 measures the data of the ocular optical 
system, such as ocular aberrations or Zernike coefficients, on 
5 the basis of the measured received light signal (4) , (7) and/or 
(10) (S105). The arithmetic part -210 performs a correction 
image simulation (S107). The details will be described later. 
The arithmetic part 210 outputs data to the correction image 
simulation data display part 230 and the memory 24.0 (S109). 
10 Fig. 5 is a flowchart of the correction image simulation 

of the step S107. 

The arithmetic part 210 calculates a best image condition 
(S201). As described later, the details are such that the 
arithmetic part 210 obtains a lower order Zernike coefficient 
15 so that the Strehl ratio becomes maximum or the phase shift 
becomes as small as possible, and obtains corrective correction 
data. As the corrective correction data, suitable data can be 
named among, for example, coefficients corresponding to 
defocus, astigmatism components, S, C, A, higher order 
20 spherical aberrations, higher order astigmatism aberrations, 
higher order coma aberrations, the Strehl ratio and the like. 

The arithmetic part 210 obtains the wavefront 
aberrations W(x, y) at the time of the best image condition, 
and calculates the pupil function f (x, y) from W(x, y) by the 
25 following expression (S203) . 



20 
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W(X,Y) = Y£c*>->Z>'-<{X,Y) (111. iSiSi) 

f (x, y) =e ikW(x ' *>' 

The arithmetic part 210 calculates the luminous 
distribution function Land(x, y) of the Landolt's ring (or 
arbitrary image) with reference to the memory 240 (S205). The 
5 arithmetic part 210 performs a two-dimensional Fourier 
transformation to obtain a spatial frequency distribution FR(u, 
v) (S207). The arithmetic part 210 obtains a frequency 
distribution OR(u, v) after passing through the ocular optical 
system by multiplying the spatial frequency distribution FR(u, 
10 v) of the Landblt's ring (or arbitrary image) and an ocular 
spatial frequency distribution OTF(u, v) together as indicated 
by the following expression (S209) . 

FR(u, v) x OTF(u, V) -> OR(u, v) 

Next, the arithmetic part 210 performs a two-dimensional 
15 inverse Fourier transform to obtain a luminous distribution 
image Landlmage (X, Y) of the Landolt's ring (or arbitrary 
image) (S211) . The arithmetic part 210 displays the LandImage(X, 
Y) and PSF(X, Y) on the display part 230 by a suitable display 
method of a drawing, graphic data, a graph and/or a numerical 
20 value, and suitably stores the data in the memory 240 (S213). 
The arithmetic part 210 reads but corrective correction data 
from the memory 240 as the need arises, and outputs it to the 
display part 230 (S215). 

25 Hereinafter, a first and a second examples of a detailed 
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flowchart concerning the step S201 will be described. 

Fig. 6 shows the flowchart concerning the first example 
of the best image condition calculation. 

- First, the arithmetic part 2io sets a threshold value for 
5 respective aberration quantities RMSi 23 " 1 as a branch condition 
(S401). For example, this threshold value can be made a 
sufficiently small value (for example, 0.1) of aberration. The w 
arithmetic part 210 calculates the Zernike coefficients Ci 23 " 1 
from the measured detection wavef ront , and converts them to the 
10 aberration quantities RMSi 2 *" 1 by the. following expression 
(S403). 

I 2j-i 

RMSf J -' = , C? J -' 

(s? J -' = 2 (2 J = 0 , e? J -' = 1 (2j * 0 ) 

The arithmetic part 210 Judges whether at least one of 
the values of RMSi 2 *" 1 (i > 2) is the threshold value or higher 
15 (S405). Here, in the case where a Judgment of No is made, it 
proceeds to step S419. On the other hand, here, when a judgment 
of Yes is made, the arithmetic part 210 carries out a next 
processing. 

That is, the arithmetic part 210 judges whether at least 
20 one of the higher order spherical aberration quantities R 4 °, 
R 6 ° ••• is the threshold value or higher (S407) . Here, in the case 
of Yes, the arithmetic part 210 causes the aberration to change 
a coefficient (C 2 °) corresponding to the defocus so that the 
Strehl ratio becomes maximum (S409) , and on the other hand, in 
25 the case of No, it proceeds to step S411. Next, the arithmetic 
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part 210 judges whether at least one of the asymmetrical higher 
order coma aberration quantities RMSi 2 - 1 " 1 (i: odd number) is the 
threshold va^e or higher (S411) • H&pe, in the case of Yes, the 
arithmetic part 210 causes the aberration to change the 
5 coefficient (C 2 °) corresponding to the defocus so that the 
Strehl ratio becomes maximum (S413) , and on the other hand, in 
the case of No, it proceeds to step S415. Further, the 
arithmetic part 210 judges whether at least one of the higher 
order astigmatism aberration quantities RMSi 2 - 3 " 1 (i: even number 
10 and 2j-l * 0) is the threshold value or higher (S415). Here, 
in the case of Yes, the arithmetic part 210 adds astigmatism 
components (C 2 ~ 2 , C 2 2 ) to the aberration so that the Strehl ratio 
becomes maximum (S417), and on the other hand, in the case of 
No, it proceeds to step S419. 
15 in this way, the arithmetic part 210 calculates OTF(u, 

v) and PSF(X, Y) from the aberrations, and further calculates 
the corrective correction data (suitable data such as 
coefficients corresponding to the defocus, astigmatism 
components, S, C, A, higher order spherical aberrations , higher 
20 order astigmatism aberrations, higher order coma aberrations, 
and Strehl ratio) from the Zernike coefficients, and stores 
them in the memory 240 (S419). 

Incidentally, in order to correct only a desired 
component among the defocus and the astigmatism components, any 
25 of the pairs of the steps S407 and S409, the steps S411 and S413, 
and the steps of S415 and S417 may be omitted, or a step may 
be added to correct suitable higher order aberrations or 
Zernike coefficients other than these. For example, in the case 
where a fourth-order spherical aberration is mainly included 
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In the higher order aberrations , the corrective correction data 
can be obtained by correcting in the direction in which the 
defocus quantity corresponding to the lower order aberrations 
are increased/ 

5 Next, the detailed processing of the steps S409, S413 and 

S417 will be described. In the respective steps, the arithmetic 
part 210 carries out the processing as follows. 

The wavefront at the time of the objective complete 
correction calculated from the Zernike coefficients is 
10 expressed by the following expression. 



n i 



/=o y=o 

In order to obtain a more suitable image plane, the 
arithmetic part 210 adds to the wavefront aberrations W(x, y) 
the lower order Zernike coefficients Ci^^d £ i s; 2) at each 

15 step presently noted for the aberration quantities comparable 
to the higher order aberration quantities according to the 
threshold value of the last noted higher order aberration 
quantities (RMS 4 °, RMS 6 ° - ) in the flow. For example, C 2 ° is 
added at the step S409; C 2 °, at the step 413; and C 2 ~ 2 , C 2 2 , at 

20 the step S417. 

W{XJ)^Y^C) H Z] H (XJ) (i^2, 2SjSi) 
M> y=0 

Further, the pupil function f (x, y) is obtained from the 
wavefront aberrations in the manner described below, 
f (x, y) = e"* 1 "' y) 
25 (i: imaginary number, k: wave number vector (2k/X) , X: 

24 



WO 2004/019771 



PCT/JP2003/010690 



wavelength) 

The arithmetic part 210 performs the Fourier 
^ transformation on this pupil function f(x, y) , so that an" 
Amplitude distribution U(u, v) of a point image is obtained as 
5 in the following expression. 



00 

amplitude 

U(u,v) = J Jf (x,y) exp 



i 2^/ \ 
(ux + vy) 



R X 



dxdy 



(X: wavelength 

R: a distance from a. pupil to an image point (retina) 
(u, v) : a coordinate value on a plane orthogonal to an optical 
10 axis while an image point O is made the origin 
(x, y) : a coordinate value on a pupil plane) 

The arithmetic part 210 multiplies U(u, v) by its complex 
conjugate and obtains I(u, v) as a point image intensity 
distribution (PSF) by the following expression. 
15 I(u, v) = U(u, v)U*(u, v) 

Besides , when the center Intensity of PSF at the time when 
there is no aberrations (W(x, y) = 0) is made Io(0, 0), the 
Strehl ratio is defined as follows: 

Strehl ratio - 1(0, 0)/Io(0, 0). 
20 In the first example, the arithmetic part 210 recursively 

or analytically obtains a value of the lower order Zernike 
coefficient C 13 (1 =s i s 2) so that the value of the Strehl ratio 
becomes maximum. 

Next , Fig . 7 is a flowchart concerning the second example 
25 of the best image condition calculation. 

First, the arithmetic part 210 sets a threshold value for 
the respective aberration quantities RMSi 23 " 1 as a branch 
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condition (S501) . For example, this threshold value is made a 
sufficiently small value (for example, 0.1) of aberration. . 

The arithmetic part 210 calculates the Zernike 
coefficients Ci 2j-i from the measured detection wavefront, and 
5 converts them to the aberration quantities RMSi 23 ' 1 by the 
expression indicated' in the first example (S503). "The 
arithmetic part 210 judges whether at least one of RMS i 23 " 1 (i 
> 2) is the threshold value or higher (S505) . £ere; in the case 
of the judgment of No, it proceeds to step S519. On the other 

10 hand; in the case of the judgment of Yes, the arithmetic part 
210 carries out a next processing. 

That is, the arithmetic part 210 judges whether at least 
one of the higher order spherical aberration quantities R 4 °, 
R 6 ° ••■ is the threshold value or higher (S507) . Here, in the case 

15 of Yes, the arithmetic part 210 causes the aberration to change 
the coefficient (C 2 °) corresponding to the defocus so that the 
phase shift becomes as small as possible (S509), and on the 
other hand, in the case of No, it proceeds to step S511. Next, 
the arithmetic part 210 judges whether at least one of the 

20 higher order coma aberration quantities RMSj 23 " 1 (i: odd number) 
is the threshold value or higher (S511) . Here , in the case of 
Yes, the arithmetic part 210 causes the aberration to change 
the coefficient (C 2 °) corresponding to the defocus so that the 
phase shift becomes as small as possible (S513), and on the 

25 other hand, in the case of No, it proceeds to step S515. Further, 
the arithmetic part 210 judges whether at least one of the 
higher order astigmatism aberration quantities RMSi 2 ^ 1 (i: even 
number and j * 0) is the threshold value or higher (S515) . Here, 
in the case of Yes, the arithmetic part 210 adds the astigmatism 
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components (C 2 "* 2 , C 2 2 ) to the aberration so that the Strehl ratio 
becomes maximum (S517), and on the other hand, in the case of 
No, it proceeds to* step S519. 

In this way: the arithmetic part 210 calculates OTF(u, 
5 v) and PSF(X, Y) from the aberrations, and further calculates 
the corrective correction data (suitable data such as the 
coefficients corresponding to the defocus, astigmatism 
components, S, C, A, higher order spherical aberrations , higher 
order astigmatism aberrations , higher order coma aberrations , 
10 and Strehl ratio) from the Zernike coefficients, and stores 
them in the memory 240 (S519). 

Incidentally, any of the pairs of the steps S507 and S509, 
the steps S511 and S513, and the steps S515 and S517.may be 
omitted so that only a desired component is corrected among the 
15 defocus and the astigmatism components . Besides , a step may be 
added so that suitable higher order aberrations or Zernike 
coefficients other than these is corrected. 

Next, the detailed processing of the steps S509, S513 and 
S517 will be described. The arithmetic part 210 carries out the 
20 processing as follows. 

First, as described in the detailed processing of the 
steps S409, S413 and S417, the arithmetic part 210 obtains the 
point image intensity distribution (PSF) from the expression 
of the wavefront at the time of the objective complete 
25 correction calculated from the Zernike coefficients. Next, the 
arithmetic part 210 performs a Fourier transformation (or 
autocorrelation) on the PSF to normalize it as in the following 
expression and obtains OTF. 
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00 



R(r,s) = J Jl(u,v)e- i2 ^ (ru+8v) dudv 



—00 



(r, s: a variable of a spatial frequency region) 




In general, the amplitude of a spatial frequency region 
and a phase distribution R(r, s) become complex numbers, and 
when its real number part is A(r, s), and its imaginary part 
5 is B(r, s) , 

R(r, s) = A(r, s) + iB(r, s) 

and the shift of the phase (phase shift, PTF) is 
calculated by 



recursively and analytically obtains such a value of the lower 
order Zernike coefficient d 23 " 1 that a value at which the R(r, 
s) has an extreme value is brought to a high frequency to the 
extent possible, that is, the phase shift becomes as small as 
15 possible. 

Incidentally, with respect to the first example and the 
second example of the best image condition calculation, both 
the processings may be carried out to obtain such a condition 
that the Strehl ratio is large, and the phase shift is small. 

20 

5. Display example 



(f>(r\s) = tan 



10 



In the second example, the arithmetic part 210 
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With respect to the best image display - Strehl ratio 
optimization (pupil diameter of 4 mm). Fig. 9 displays , as 
numerical data, sphere S, cylinder* C, astigmatism axis angle 
Ax of the corrective correction data, and sphere S, cylinder 
5 C, and astigmatism axis angle Ax of measurement values before 
corrective correction. In this example, since higher order 
aberration components have specified values or higher, there 
occurs a difference between the corrective correction data and 
the measurement values. 
10 With respect to the best image display - PTF optimization 

(pupil diameter of 4 mm). Fig. 9 displays, as numerical data, 
sphere S, cylinder C, and astigmatism axis angle Ax of 
corrective correction data, and sphere S, cylinder C, and 
astigmatism axis angle Ax of measurement values before^ 
15 corrective correction. In this example, since a higher order 
aberration components have specified values or higher, there 
- occurs a difference between the corrective correction data and 
the measurement values. 

In these drawings, the, wavef ront aberrations, PSF, OTF, 
20 OTF (two-dimensional display), S, C, Ax, Landolt's ring, 
appearance of index, and the like are displayed on the display 
part. 0 

Fig. 10 is a view showing a display example concerning 
the comparison (pupil diameter of 4 mm) of pre -correct ion and 
25 post-correction. In this drawing, the wavefront, the 
appearance of the Landolt's ring, and the Strehl ratio before 
the correction and those after the correction are displayed. 
As shown in the drawing, it is indicated that after the 
correction, the Strehl ratio is higher, the wavefront 
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aberrations become relatively uniform, and the Landolt 1 s ring 
is also relatively well seen. 

Fig. 11 is an explanatory view of an example of 
prescription data (pupil diameter of 4 mm) for eyeglasses or 
5 contact lenses. Fig. 12 is an explanatory view of an example 
of data (pupil diameter o£ 4 mm) for a refraction correcting 
surgical operation. 

These data are stored in the memory 240 by the arithmetic 
part 210, and/or displayed on the display part 230. This example 
10 indicates that in the data of a case where the refraction 
correcting surgical operation is performed while only SCA are 
made correction data, by performing correction to intensify the 
value of S in the correction data, to weaken the value of C, 
and to slightly change the axial direction of A, the Strehl 
15 ratio becomes high and the correction effect becomes high. 

A modified example of the invention will be described 

below. 

This modified example modifies the calculation method of 
the best image condition at S201 of Fig. 5. 
20 A component of an i-th row and a j-th column of Jacobian 

matrix A is 

A _ df, 0 ) 

Where, f±(x) Is the Strehl ratio, the PTF corresponding 
to a suitable frequency, or some values of the PTF corresponding 
25 to plural frequencies. Besides, it may be a combination of the 
Strehl ratio and the PTF. Besides, a vector x is an adjustable 
parameter, and here, the sphere (or defocus corresponding to 
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that) and two astigmatisms correspond to that* 

The calculation expressions of the Strehl ratio and the 
PTF are already given . The ideal value of the Strehl ratio is 
1. It is assumed that the following expression expresses the 
5 Strehl ratio, 

,fi(x) = fs(hc, c 2 °, c 2 ~Vc 2 2 ) (2) 
Where, fi denotes the expression of the same indication 
in the^expression (1). 

Besides, for example, as the PTF, values corresponding 
10 to the spatial frequency of 3 cpd, 6cpd f 12 cpd, and 18 cpd are 
taken, and it is ideal that this is 0, 

f 2 (X> = fpTF3(hC, c 2 °, c 2 " 2 , c 2 2 ) (3) 
f 3 (x) = fpTF6(hC, c 2 °, c 2 " 2 , c 2 2 ) (4) 
f 4 (x) = f P TFi2(hc, c 2 °, c 2 - 2 , c 2 2 ) (5) ^ 

15 f 5 (x) = fpTF18(llC, c 2 °, c 2 " 2 , c 2 2 ) (6) 

In the expressions (2) , (3) , (4) , (5) and (6) , he denotes 
a vector of higher order aberration coefficients, c 2 ° denotes 
a coefficient of a def ocus term relating to the sphere, c 2 " 2 and 
c 2 2 denote coefficients of terms relating to astigmatism. The 
20 vector he is given by wavefront aberration measurement, and 
here, it is constant. Thus, the remaining three coefficients 
are made a parameter vector x and are suitably moved to guide 
f s to the minimum value, which is a task here. 

Here, the partial differentiation of the expression 1 can 
25 be calculated by slightly moving the parameters to prepare a 
change table, and the Jacobian matrix in this system is 
obtained. 

Now, when the task here is expressed in other words again, 
since nonlinear optimization in the case where the Jacobian, 
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that is , the partial differential coefficient is known has only 
to be performed, when optimizing algorism of a Newton method 
system is used, it is easy to obtain a solution since the example 
is simple. When a specific solution according to a corrected 
5 Marquardt method is stated, a correction vector Ax can be 
obtained by 

(A*WA + XI) Ax - A*W(y - f(x)) (7) 

Here, t at the shoulder of the matrix denotes a transposed 
matrix, and W denotes a weighting matrix. The first element of 
10 y corresponds to Strehl ratio, and the remainder corresponds 
to four components of PTF , it has only to be made (1,0,0,0,0)*. 
K is called a damping factor, and it is made large at first, 
and then, it is made small in accordance with going of 
optimization. . 
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(8) 



A subscript corresponds to a subscript of f . Weighting 
suitable for the object of a prescription can be freely 
performed, for example, when the Strehl ratio is desired to be 
selectively optimized, Wi is made large. The expression (7) is 
20 applied several times, and when 

S = W(y - f(x)) (9) 

becomes suitably small (when a conversion condition is 
satisfied) , the calculation is stopped, and x at that time is 
made the solution. By this, the optimum the sphere (or defocus 
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corresponding to that) and two astigmatisms are obtained. 

Incidentally, the best image condition calculation can 
also be performed by finding out the position where the 3trehl 
ratio becomes maximum or the phase shift (PTF) becomes 
5 substantially zero while the defocus amount and/or the 
astigmatism component is changed slightly. Besides, a position 
where the Strehl ratio becomes maximum or the phase shift (PTF) 
becomes substantially zero may be obtained by using a 
well-known Newton method. 

10 

6. Addition 

The correction data measurement method of the invention 
or the apparatus and system of the correction data measurement 

15 can be provided by a correction data measurement program for 
causing a computer to execute the respective procedures , a 
computer readable recording medium recording the correction 
data measurement program, a program product including the 
correction data measurement program and capable of being loaded 

20 in an internal memory of a computer, a computer, such as a server, 
including the program, and the like. 

Although the measurement data indicating the refractive 
power distribution of the subject eye is obtained by the optical 
system 100 shown in Fig. 1, the invention is not limited to this, 

25 but may be structured by another aberrometer or the like. 

Industrial Applicability 

As described above, according to the invention, from the 
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result measured by the eye characteristic measuring apparatus 
capable of performing the measurement up to higher order 
aberrations, optical performance is evaluated by, for example, 
the Strehl ratio or the phase shift in not only the case of only 
5 the higher order aberrations corresponding to the objective 
complete correction, but also the case in which the lower order 
aberrations are added, such lower order aberration quantities 
that for example, the Strehl ratio is large and/or the phase 
shift becomes small is calculated, and the correction data of 
10 S, C, A and the like at that time are obtained, so that the resuit 
closer to the subjective value can be obtained. 
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